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The comments raised in Ref. [J] by Mishra et al aim at two papers contained in Ref. 0]. We 
show that those comments on Ref. Q] pointed out by Mishra et al in Ref. [l[ are not relevant and 
the concept used in Ref. 0] is consistent and in compliance with the classical approximation of the 
transport coefficients 0]. We would also like to note that most of the comments in Ref. [J] were 
meant for light quarks, but are not even appropriate for heavy quarks. 
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The comments of Mishra et al in Ref. [J| arc organised as a list of observations ( "lacunas" ) . Here we reply to these 
observations one by one. 

Lacuna 1: Notations 

The notations and/or expressions can easily be understood depending upon the physical considerations (massless 
or massive quarks) they are applied to. In the case of high momenta (at least p > 5 GcV/c), there is, in fact, no 
significant difference between E and p, for a charm quark jet. 



Lacuna 2: Fokker-Planck (FP) Equation 

We begin with the one-dimensional FP equation 0, Q as 



^D{p,t) = ^^[T,{p)D{p,t)] + ^[r2(p)D(p,t)], (1) 

that describes the evolution of the momentum distribution, D(j),t), in a domain (— cxd < p < oo), at a given time 
t (0 < t < oo), of a test particle undergoing Brownian motion. ri(p) and T2{p) are known as the moments or the FP 
coefficients or the transport coefficients. Usually ^i{p) is related to the collisions whereas T2{p) to the momentum 
diffusion in the medium when a test particle undergoes Brownian motion. 

We now outline the classical approximations [1, Q for the transport coefficients in which the drag force. A, is 
assumed to be related to the collisional energy loss, —dE/dL, as 

dE 

r^{p)=^~—^pA{p)^pA, (2) 

where A = {A{p)) = (^^^): ^i^d the diffusion coefficient is related to the drag as 

dE 

r2(p) « pAip)p « TAip)p « T{- — } EE Vf. (3) 

Within this approximation the drag, A, and the diffusion. Dp, coefficients are momentum independent. Now one 
can write the FP equation in (P) as Q 

= A^[pDip,t)] + VF^D{p,t). (4) 
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The authors of Ref. [l| also agree to the above form (see, e.g, eq. (C3) of Ref. [l|), which uiiphes that they, obviously, 
accept the momentum independence approximation of A and Dp. 

We now quote the solution of the FP equation ^ for a given time t as obtained in Ref. 0] 



Dip,t) 



exp 



Wit) 



(5) 



where W{t) is given by 



W{t)= Uj^VF{t')exp 2 J A{t")dt" dt'j exp (^-2 A{t') dt 



(6) 



We further clarify that for a given initial condition, S{p — po), and for the momentum independence approximation^ 
of A and Vp the above solution is correct and unique which could easily be checked as the detailed calculational steps 
arc given in Ref. . It is also worthwhile to note that for massless quarks the analytically obtained solution of the FP 
equation in ([4]) agrees well with the solutions [5j if one solves ([T]) numerically for a given temperature (T = 400 McV) 
and initial momentum (po = 16 GeV/c), up to i ~ 8 fm/c, beyond which the classical approximations for transport 
coefficients are not valid. This time regime is quite appropriate to study the energy loss probability distribution of a 
jet in quark-gluon plasma (QGP) whose life time is of the same order for RHIC energies. 

Now, the authors of Ref. [l| claim that there should be also a linear term in momentum, {p — (p)), in the exponential 
of ([5]) for reducing to the massless Maxwellian form in the thermodynamic limit. Within the given initial condition, 
S{p ~ po), and the momentum independence approximation of A and "Dp, such a linear term in the exponent of 
([5|) does not appear (see also Lacuna 4). Of course, if one changes the initial condition or relaxes the momentum 
independence of A and T>p, one will then arrive at an altogether different solution. Note that the authors of Ref. [H 
have not defined clearly what they really mean by "the possibility of a different type of solution" . However, we intend 
to discuss various possibilities later in lacuna 4 while analysing the asymptotic forms of the FP equation in detail. 

Lacuna 3: Probability function in p vs. E 



According to ([5|) the function D{p, t) is a probability distribution in p in one dimension (— oo < p < oo) for a given 
time t. Since it is an even function in p the normalisation condition can be written as 



/oo /'OO pOO 

D{p,t)dp = 2 D{p,t)dp^2 f{E,t)dE^l, 
-oo Jo Jo 



(7) 



with E = p = \p\. Thus, the normalisation requirement, for the massless case, is unambiguously preserved. This also 
holds for charm quarks as long as we consider momenta p much larger than the charm quark mass. 

Lacuna 4-' Unphysical asymptotic form 



Let us make another 'thought situation' where the drag (A) and the diffusion {T)p), are time independent, i.e., the 
test particle will experience a constant amount of drag and diffusion over a given length of medium or time traversed. 
This consideration is, however, without any loss of generality for the case of a static QGP. 

In such a static medium the distribution in ^ reduces to 



Dip,t) 



A 



2ttVf (1 - 



^(1 



Now in the thermodynamic limit {t — !■ oo), the above equation simply becomes 

D{p)-- 



A 



2ttVp 



(8) 



(9) 



^ It is interesting to note that two of the present authors in Ref. 01 also used the same approximation and the said solution for a chemically 
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At this point, wc must check the uniqueness of the solution in For the purpose, we refer to the second pair of 
the characteristic equation in Ref. 0|(See. as for example equation (22) in the first article in Ref. 0) corresponding 
to ^ in Fourier space, which amounts to the i — > oo limit, that reads 

dx db 

^ ^ - (10) 
Ax Vpx^D 

where D = D{x) is the Fourier transform of the momentum distribution function. It is easy to show that the solution 
obtained using pO)l agrees perfectly with that of given in ([9]). So, the solution obtained in Ref. is unique and a 
linear momentum term in the exponential does not arise, once the FP equation is given by (C3) of Ref. Q or by 
However, we are aware of the fact that relaxing our momentum independence approximation of A and 'Dp, will 
certainly lead to different solutions, as the form of the FP equation will then be different from ^ or (C3) in Ref. 
We analyse various possibilities in the Appendix. 

Based on our analysis in the Appendix we note that the FP coefficients, Fi and F2, cannot simply be assumed in 
any form to ensure the correct form of the (equilibrium) distribution because this requires at least an exact evaluation 
of the drag (A) and the diffusion {"Dp) coefficients, which is even indeed a very involved task^ Q. But knowing 
the collisional energy loss in terms of elementary collision reaction amplitudes, one can approximate the drag as in 
(PI and the diffusion as in The form of the approximation determines, a priori, the shape of the equilibrium 
distribution as discussed in the Appendix. Obviously, the momentum independence of A and Vp produces at least 
a consistent solution over the momentum range, —00 < p < 00, in contrast to the other possibilities (see Appendix). 
Again, we would also like to note that for small t the change in diffusion"^ will affect only the width and the height of 
the distribution, but not the peak position which is determined by the drag. Since we really do not wish to obtain an 
equilibrium distribution for a jet but intend to obtain an energy-lossprobability distribution over a finite length (~ 
8 to 10 fm) of the medium, ([5]) is a reasonably good approximation [3|. In addition, it was independently verified [5| 
numerically without any approximation as given in ([T]), where Fi and F2 are obtained using kinetic theory calculations. 
Nevertheless, it would have been easier for us to discuss this problem if the authors of Ref. [ij would have been more 
specific in their statement about "a different type of solution" . 

Now, the purpose of the study in Ref. [1] was the following: before wc took up these studies [1] , the dominance 
of the radiative energy loss for the phenomenon of jet quenching in heavy-ion collision seemed to be well established 
in the heavy-ion community. Within this simple approach we reconsidered, in 2003, the role of collisional partonic 
quenching and showed that it could be significant and cannot just be overlooked as it were done in the literature. 
However, only recently, after publication of new data on the nuclear suppression factor in RHIC, this simple idea has 
gained wider interest. In fact, an additional contribution to the partonic energy loss appears to be necessary and a 
collisional component is a welcomed remedy, as advocated in 0]. At this point, since RHIC BNL has provided very 
accurate data and LHC CERN will be operational soon, one indeed needs to improve our simple approach in different 
possible ways, which would definitely be a very desirable to verify the importance of the collisional energy loss. 



Lacuna 5: Mean energy and loss 



According to our explanations on the above three points (Lacuna 2-4) the calculation of the mean energy obtained 
in Ref. Q follows from the expression of the distribution function itself in the classical approximation. The mean 
momentum of the test particle, as for example we refer to ([5]), is (p) = poe~^*, along with the diffusion process in 
momentum space as (p^) — (p)^ — 2^(1 — e~^-^*). 

However, just to be sure, one can numerically compare (p) = poe with the one given in (Cll) of Ref. They 
agree up to three to five decimal places for time interval (1 — 10) fm/c as obviously the second term in (Cll) is 
negligible and our classical approximation is consistent. This suggests that using (p) = poe"'^*, is not a matter of 
great concern keeping in mind again the life-time of quark-gluon plasma. However, we note that it was used only for 
light quarks but proper care has been taken for heavy quarks. 



Lacuna 6: Convolution integral for hadron spectrum 



We could not make out the point raised by the authors of Ref [l| on our numerical calculations. In particular 
when they talk about a fixed length, L, in the convolution. However, the length L{(j)) as defined in Ref. [3| will be 



^ Here, we do not even wish to speak of the higher order FP eoefHcients in the Taylor expansion. 
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determined by a jet created at the transverse position, r and the production angle, cj) in central collision. With this L 
one should convolute the spectra along with the geometry given by a cylinder of radius R in Bjorken hydrodynamics, 
which, however, should be restricted up to the critical temperature, Tc ^ 0.2 GeV. We checked our numerical results 
and reproduced Q{pt) for light hadrons. The results obtained for the scaled energy loss, for heavy quarks have 
been independently verified in Ref. Q. Hence, we believe that the approach as well as the analytical and numerical 
solutions in Ref.Q are correct as long as they are used within the domain of (8 — 10) fm/c. 
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APPENDIX A: ASYMPTOTIC SOLUTION OF FOKKER-PLANCK EQUATION WITH VARIOUS 

POSSIBILITIES OF ITS COEFFICIENTS: 

For convenience we write down the general form of the FP equation from ([T]) for steady state as 

§-^D{p,t) = ^[T,ip)D{p,t)] + ^[r,ip)Dip,t)] = 0. (Al) 

where D{p, t) is, in general, the probability distribution in momentum, p (— oo < p < oo) at a given time t (0 < t < oo) 
of a test particle. Now we investigate the asymptotic solutions of the FP equation with various choices of its coefficients, 
Ti{p) in (HI) and r2(p) in ([3]). 

Type-I: 

We consider Ti{p) « pA and r2(p) ~ Ap^. Now it is easy to show that the steady state solution of the FP equation 
in (|Aip reads as 

Dip) - ^ , (A2) 

where C'l is a constant. This solution is not Maxwellian and diverges in the limit, p 0. Thus, it is not a consistent 
solution over — cx) < p < oo. 

Type-II: 

We now consider Fi {p) « pA and T2{p) ~ ATp. Similarly, the form of the steady state solution of the FP equation 
in (jAip can be obtained as 

D{p) = ^ e-* . (A3) 
P 

This appears like a pseudo-Maxwellian form but is again not a consistent solution over — cxd < p < oo as it diverges 
in the limit, p — > 0, like the Type-I above. However, it is very easy to anticipate the reason. 

Type-Ill (our case): 

We consider Fi(p) « pA and T2 = Vp- The steady state solution is 

D{p) = e"^ , (A4) 
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Type-IV : 

We now explore the simplest one where both Fi and r2 are momentum independent. Then the steady state solution 
reads 

D{p) = d e~^P . (A5) 
Now one can claim it to be of the Maxwcllian form just by choosing = i. 
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